Collective response of self-organised clusters of mechanosensitive channels 



Ksenia Guseva/'0 Marco Thiel,^ Ian Booth, ^ Samantha Miller,^ Celso Grebogi/ and Alessandro de Moura^ 

^Institute for Complex Systems and Mathematical Biology, King's College, 
University of Aberdeen, Aberdeen AB243UE, United Kingdom 
^Institute of Medical Sciences, University of Aberdeen, Aberdeen AB252ZD, United Kingdom 

(Dated: January 4, 2011) 

Mechanosensitive channels are ion channels activated by membrane tension. We investigate the 
influence of bacterial mechanosensitive channels spatial distribution on activation (gating). Based 
on elastic short-range interactions we map this physical process onto an Ising-like model, which 
enables us to predict the clustering of channels and the effects of clustering on their gating. We 
conclude that the aggregation of channels and the consequent interactions among them leads to a 
global cooperative gating behaviour with potentially dramatic consequences for the cell. 



Recent advances in the understanding of the functional 
organisation of the cell membrane are shedding light 
on the complex dynamics of the components of the cell 
surface [1-3]. There is considerable evidence that non- 
specific, membrane- mediated forces are important for the 
formation of protein complexes on the membrane [4]. In 
this work, we address the question of how the spatial 
organization of membrane proteins can be shaped by 
their short-range interactions, and how this spatial or- 
ganization can affect their function. We focus on the 
behaviour of mechanosensitive channels, which are ac- 
tivated by membrane deformation. These channels are 
present in several organisms, such as bacteria, humans 
and plants, and are responsible for a variety of functions, 
ranging from volume regulation, locomotion, and sensory 
input and signalling [5]. Membrane mechanical proper- 
ties can also influence other types of ion channels [6] ; thus 
mechanosensitive channels can be used as a general model 
for proteins actuated by membrane- mediated forces. The 
channels which are best characterized are from E. coli, 
and are of two types: channels of large (MscL) and small 
(MscS) conductance [7]. These channels are responsible 
for preventing the osmotic pressure from reaching danger- 
ous levels under hypoosmotic shock. They are activated 
directly by membrane tension, which causes a membrane 
deformation in the channel neighbourhood [8], and gate 
(open) when the cell is placed in an environment of low 
osmolarity. Although there are not many mechanosen- 
sitive channels on the E. coli membrane, they are often 
overexpressed for studies, and are present in larger num- 
bers in other organisms j5|. For these reasons, they are 
an excellent and well-studied model system, from which 
broader conclusions may be drawn. Since the forces in- 
fluencing them are short-range, and typically attractive, 
there is the possibility that they form clusters. Two ques- 
tions of crucial importance are: (i) under what conditions 
do channels cluster? (ii) and how does clustering affect 
the gating of the channels? We address both questions 
by presenting a general statistical mechanics framework 
— which can be easily carried over to other types of chan- 
nels with different interaction forces — and conclude that 
channels should indeed cluster if their density (or interac- 



tion strength) is high enough, and that this has an enor- 
mous effect on the gating response of the channels, which 
display a much richer cooperative behaviour than can be 
inferred from the characteristics of individual channels. 

The collective behavior of MscL is a result of their mu- 
tual interactions, mediated by the membrane deforma- 
tion around them. We model channel agglomeration and 
its effects on gating in a coarse-grained manner, in which 
the channels are placed on a two-dimensional lattice and 
only their nearest-neighbour interactions are considered. 
Initially, using a lattice gas model, we obtain the condi- 
tions for channel agglomeration, and their detailed spa- 
tial configuration. Based on that, we subsequently show 
that the gating dynamics of the clustered channels can 
be mapped onto an Ising-like model, with the addition 
of a spatially inhomogeneous field. This opens a new ap- 
proach in the analysis of propagation of conformational 
states through a cluster of proteins. One of our ma- 
jor findings is that clustering leads to a lower threshold 
of channel activation, causing the clustered channels to 
open for lower membrane tensions than in the case of 
isolated channels. Furthermore, our method allows us to 
study non-equilibrium properties of the system such as 
the dynamics of transition. Due to this transition, clus- 
tering leads to an increase in the time it takes for the 
clustered channels to open in response to osmotic shock. 
Both these results show that the channel response to os- 
motic stress is crucially affected by interactions inside 
clusters. 

Membrane proteins diffuse in the lipid bilayer, which 
can be considered as a two dimensional fluid. They also 
interact with each other [1]. The possible forces among 
membrane proteins are electrostatic and membrane- 
mediated interactions j9|. The electrostatic forces can be 
neglected due to the charge screening effect in physiolog- 
ical solutions, since the Debye length is 1 nm |2(, which 
is just a fraction of the size of a channel molecule. Segre- 
gation by lipid affinity is also not considered since MscL 
does not exhibit strong lipid preferences in E. coli flQ| . 
Therefore, we consider only elastic forces in what fol- 
lows These forces are short range and have magni- 
tude of the order of k^To (Tq is the typical room tem- 
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perature, Tq = 300K). They arise from the hydrophobic 
misraatch between the size of the core of a protein and 
the length of the hpid layer that surrounds it Since 
lipids are more flexible than proteins, they tend to deform 
and adapt to the size of the protein core (see Fig. [1^). 
It is the energy cost of deforming the lipid layers which 
results in a force among nearby channels on the mem- 
brane. As shown in Fig. [1^, this kind of deformation 
depends on whether the channel is open or closed; there- 
fore, the forces between channels depend on their state. 
Since E. coli is the paradigm for such studies, we use the 
appropriate parameters for the MscL of these bacteria to 
obtain the deformation profile (/)(r), defined as the dis- 
tance of the membrane contour to its relaxed state, as 
a function of the distance r to the protein. For that we 
minimise the free energy. 
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where / is the mismatch length, the first term repre- 
sents the energy cost of membrane stretching {Ka = 
GO/c^To/nm^), second of membrane bending (/^^ = 
2O/C5T0, with Co = 0.009nm~^ as membrane curvature, 
we considered values in the range 0-0.04 nm~^, without 
any significant difference) (sl, HH, and last term con- 
siders the membrane tension r. To find the (j){r) which 
minimises G, we numerically solve the Euler-Lagrange 
equations corresponding to Eq. ([T]) (see for details). 
Considering a system of two channels, we find that the 
forces between them are short range (~ 5 nm) (Fig.flJ)), 
and they remain roughly the same for slight differences 
in value of cq and r. This distance is comparable to the 
diameter of a single protein, which validates our coarse- 
grained 2D lattice approximation, which will follow. 
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FIG. 1. (Color online) a) Examples of the deformation profile 
around two open and two closed channels respectively, b) 
Interaction energy between two channels in different states 
(as in [H).) 

The dynamics of agglomeration. We consider initially 



the normal physiological conditions for bacteria, which 
corresponds to a low membrane tension. In this case, 
all the channels are in the closed state and diffuse on 
the cell, interacting with each other through the elas- 
tic forces described above. Using the fact that their in- 
teraction is short range, and can be approximately de- 
composed into pairwise interactions, we place the chan- 
nels on a 2D square lattice, and describe the system by 
the Hamiltonian H = — ^Xl^^^^ SiSj, where Si are occu- 
pancy variables associated with the lattice sites, which 
assume values either 1 or —1 for an occupied or free lat- 
tice site respectively, and J is the strength of the pair- 
wise interaction between two channels. The brackets in 
the sum indicate that it is performed over all pairs of 
adjacent sites. Since the number of channels is con- 
stant — Yli — ^ Pi where N is the total number of 
sites, and p is the density of occupied sites — this cor- 
responds simply to a lattice gas, which is characterised 
by the existence of either a non-homogeneous or a ho- 
mogeneous phase, depending on the interaction strength 
J (or on the thermal fluctuation T), and on the particle 
density p, as follows: For a weak interaction (or a high 
temperature) only a homogeneous phase is observed in- 
dependently of the density. For J high enough (or T 
low enough) this system has two critical densities [lH, 
p± = ^(1 ± (1 - csch^(2/3J))i). When the density p is 
lower then p_ or higher then p+ the particle distribution 
is homogeneous. However, if p_ < p < p+ the parti- 
cles segregate in different domains with different local 
densities: a low density (p_) region and a high density 
(p+) region (Fig. [2]). From Eq. [TJ we can determine that 
J ^ I.25/C5T0, which gives p- = 1.710~^ particles/site. 
Given that an wild type E. coli cell has on average only 5 
channels |^ and a membrane area of -6.10-12 m^ hj, 
the density of the channels is given by p ~ 1.6 x 10"^ 
particles/site, which corresponds to a very low density, 
deep in the homogeneous phase, without any clustering. 
Below, when we consider the effect of agglomeration on 
gating, we will argue that there is a likely biological rea- 
son for this. Often, these channels are artificially over- 
expressed, and some authors do see non-homogeneous 
distributions in such situations [13, EH, suggesting that 
clustering may be occurring. Unfortunately they do not 
estimate the channel density. We suggest as a possible 
experiment to estimate the density in those samples and 
compare the appearance of agglomeration changing the 
density conditions. 

Gating dynamics. We turn now to the gating response 
of the channels, when the tension is changed abruptly. 
This is the case, for example, if the osmolarity of the 
medium is suddenly decreased. We note that the gating 
dynamics take place on a shorter time scale than the dif- 
fusion of the channels: the gating response of channels 
is of the order of microseconds [16| and the free diffusion 
is of the order of ~ {).bnm? / ps for crowded environment 
such as a biological membrane Since the area of a 
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FIG. 2. (Color online) Homogeneous and non- homogeneous 
state phase diagram, for J ■ 
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single channel is approximately 7r(2.5)^ ~ 20nm^, they 
cannot move significantly during a gating event, so we 
assume simply that their positions remain fixed in their 
initial values given by the lattice gas model discussed 
previously. Then the question is: how does the spatial 
clustering affect the channel's response to osmotic ten- 
sion? We will describe the state of each channel i by a 
variable cr^ , which can have values 1 and — 1 , correspond- 
ing to an open or closed state, respectively [18]. The new 
energy of the system can now be written as the sum of 
the non-interaction energy of the individual channels and 
their interaction energy, H = H^on + ^int • These ener- 
gies can be obtained by solving Eq. (TJ for a system of 
only one or two channels, respectively, and considering 
all the different channel conformation states. This re- 
sults in i^non = h (Ji , where 



/i = (AGgate-rAA)/2, 



(2) 



is a global non-interaction field, where AGgate ~ SO/c^Tq 
is the energy difference between conformations, A A is 
the deformation area of the protein, and r is the the 
membrane tension, which changes according to the os- 
molarity of the medium. In an analogous fashion, we 
can obtain the energy levels for interacting channels, as 
shown in Fig. [3l for an approximate distance of around 
6-7 nm from their centers. Since the interaction ener- 
gies do not change significantly with tension [l3|, we 
have assumed tq ~ 2.5ki}To/nm'\ which is the tension 
for which a single channel opens [7|. These energy levels. 



E open-closed ~ lOkbTo 
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FIG. 3. Interaction energies for two channels. For E = lOkbTo 
and P = lOkbTo 

can be written as a Hamiltonian composed of a sym- 
metric term, plus an additional spatially correlated field. 



tial correlation of the second term is due to the fact that 
it is summed only over neighboring sites. The second 
term can be interpreted as a non-interaction Hamilto- 
nian, with particles subject to a local field /c^, which is 
equal to the number of occupied neighbours of site i, 
which results in — ^ kiCii. Thus, the complete energy 
of the system is given by the Hamiltonian 
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where the sums are taken over occupied sites, and neigh- 
bouring sites. Due to the presence of the local field /c^, 
this is not simply an Ising Hamiltonian. It bears resem- 
blance to the random field Ising model [l9[, where the 
values of ki are replaced by randomly distributed local 
variables. In our case, however, they are not random, 
but instead they represent a quenched correlated disor- 
der which is a byproduct of the diffusion and aggregation 
of the channels, as modelled by the lattice gas dynamics 
discussed previously. 

The system given by the Hamiltonian in Eq. [3] is inves- 
tigated using Monte Carlo simulations. We choose bio- 
logically relevant initial conditions in which all the par- 
ticles are initialised in the closed state (mimicking con- 
ditions before an osmotic shock). This corresponds to a 
met ast able state for the interacting particles: The tran- 
sition to a global minimum (all channels open) involves 
the particles leaving a local minimum of the energy (all 
channels closed), and temporarily assuming anti-aligned 
states with respect to their neighbours. This process can 
be extremely slow, and it determines the response time 
of the bacterium to osmotic shock. This is potentially 
an important aspect of the system's dynamics: cluster- 
ing can cause the channels' response time to increase. 
We initially examine the two extreme cases where the 
particle distribution is uniform. In the low-density limit 
{p < p_), most particles are isolated and there is no ag- 
glomeration. Without any mutual iteration, the individ- 
ual channels will assume a preferred state given directly 
by Eq. [2l and will thus open when rAA > AGgate- The 
other extreme case is when the lattice is completely cov- 
ered by particles, i.e., p = 1 and ki = A for all i. In 
this case, Eq. [3] becomes the standard Ising model with 
external field h — 2P. The term 2P is due to the chan- 
nel interactions and, as a result, the transition to open 
channels occurs for a lower tension than in the case of 
non-interacting particles [P = E = {)). 

For the intermediate case, when p_ < p < the pro- 
cess of diffusion described previously makes the channels 
agglomerate in clusters of finite size with highly irregu- 
lar geometrical structures. These structures introduce an 
anisotropy in the local field /c^, which enables certain con- 
figurations of mixed states —1 and +1 to coexist. This 
can be observed by a comparison of the gating threshold 
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for a group of channels inside clusters with ramified (e.g. 
J = 0.75) or with dense (e.g. J = 1.25) structures, as 
in Fig. m The transition between the two uniform global 
states can occur in several steps, where certain discrete 
groups of channels change their configuration at different 
values of tension. The most obvious steps correspond to 
the high-density regions {ki = 4), and the low-density 
regions {ki = 0), which gate for lower and for higher ten- 
sions, respectively. However there are also intermediate 
transition steps that correspond to the outer layer of the 
clusters (particles with 1 < ki < 4), and the number 
of such steps varies with density. In the case of com- 
pact clusters the transition for ki = 4 is the only one 
present. If we consider now a system starting from the 
metastable state where all channels are closed, after a fi- 
nite (but large) simulation time of 10^ Monte-Carlo steps 
per particle, the situation changes significantly, as can be 
seen in Figs.lU^b) and Figs.lU^d). In this case, the several 
discrete transition steps are replaced by a single global 
transition at significantly larger values of tension, char- 
acterizing a delay in the reaction time of the channels. 
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The gating response time of a non-interacting channel 
is of the order of a few microseconds |l6|, a nd the sur- 
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FIG. 4. (Color online) Average conformation {a) as a function 
of particle density p and membrane tension r. The top row 
corresponds to for J = 1.25KbTo and the bottom row to J = 
0.75^6^0. (a) and (c) Equilibrium configuration; (b) and (d) 
Initial condition a = — 1 after a transient of 10^ iterations per 
particle; Simulations were carried out on a square lattice of 
linear size 1000 and averaged over 5 independent realizations. 



We now turn to a more detailed comparison to the E. 
coli system. In bacteria the physiological conditions cor- 
respond to a membrane tension of r ~ 
In this condition the independent channels are all closed, 
since their gating threshold tq is ^ 2.5k})To/nm'^ . How- 
ever, our theory predicts that in this situation r is enough 
to trigger the gating response of clustered channels (see 
Fig. [11(a)); this is a dramatic change caused by cluster- 
ing on the gating dynamics 0|. However, because of 
the initial metastable state, clustering also increases the 
delay in the gating response. If long enough, this de- 
lay may be crucial, since the bacterial cell wall can only 
sustain high tensions for a very limited amount of time. 



vival time of the order of 100 milliseconds [2l||. Using the 
approximation derived in [111, we estimate the response 

time for p = 1 to be t ~ ^k^T(2P-h) Monte-Carlo steps 
(mcs) per particle. Assuming that each mcs per particle 
correspond to the characteristic reaction time of a single 
channel, i.e., ~ 3/is, this gives us a gating response of 
^ 20 milliseconds, for r = tq. Thus we predict that un- 
der these conditions most cells would survive the osmotic 
shock even with clustered channels, but their gating time 
is orders of magnitude greater than assuming isolated 
channels. This is a measurable effect with current exper- 
imental techniques, and this delay could be crucial for 
cells with weakened cell walls, for example. Also, even 
if the channels manage to open in time, the formation 
of compact channel clusters will cause problems for the 
closing of the channels (for which there will also be a de- 
lay), after the osmotic stress is removed. Clustering has 
therefore dramatic consequences for the non-equilibrium 
dynamics of the gating. 
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